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Abstract — Nine new 2-D OOCs are presented here, all sharing 
the common feature of a code size that is much larger in relation 
to the number of time slots than those of constructions appearing 
previously in the literature. Each of these constructions is either 
optimal or asymptotically optimal with respect to either the 
original Johnson bound or else a non-binary version of the 
Johnson bound introduced in this paper. 

The first 5 codes are constructed using polynomials over finite 
fields - the first construction is optimal while the remaining 4 
are asymptotically optimal. The next two codes are constructed 
using rational functions in place of polynomials and these are 
asymptotically optimal. The last two codes, also asymptotically 
optimal, are constructed by composing two of the above codes 
with a constant weight binary code. 

Also presented, is a three-dimensional OOC that exploits the 
polarization dimension. 

Finally, phase-encoded optical CDMA is considered and con- 
struction of two efficient codes are provided. 

Index Terms — Optical orthogonal codes, optical CDMA, 
OCDMA, two-dimensional codes, 2-D OOC, wavelength-time 
hopping codes, Johnson bound, phase-encoded OCDMA. 



I. Introduction 

There has been an upsurge of interest in applying code 
division multiple access (CDMA) techniques to optical net- 
works - optical CDMA (OCDMA) [1]. This is partly due to 
the increase in security afforded by OCDMA (as measured, 
for instance, by the increased effort needed to intercept an 
OCDMA signal) and partly due to the flexibility and simplicity 
of network control afforded by OCDMA. 

There are two main approaches to data modulation and 
spreading in optical CDMA (OCDMA). The first approach, 
known as direct-sequence encoding [1], makes use of on- 
off-keying (OOK) data modulation and unipolar spreading 
sequences with good correlation properties. Traditionally, as 
in the case of wireless communication, the spreading has 
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been carried out in the time domain and we will refer to this 
class of OOCs as one-dimensional OOCs (1-D OOCs) [2]- 
[25]. A drawback of 1-D OOCs is the requirement of a large 
chip rate. By employing two-dimensional optical orthogonal 
codes (2-D OOCs) that spread in both time and wavelength 
domain, it turns out that the large-chip-rate requirement can 
be substantially reduced. There is a considerable literature on 
2-D OOC constructions. However, in this paper, we focus our 
attention only on optimal and asymptotically optimal construc- 
tions. A quick overview of optimal (or asymptotically optimal) 
constructions of 2-D codes in the literature is presented in 
Table [I] The code constructions presented here came about 
as a result of a DARPA-funded project [26], a subgoal of 
which was coming up with constructions that would give the 
experimentalist the maximum possible flexibility in choosing 
the parameters A and T. In this context, Fig. [2] provides a 
visual depiction of the parameter sets for small A, T in the 
range, 2 < A < 17, 2 < T < 33 for which constructions are 
now available as a result of the 9 constructions presented here. 

The second OCDMA approach uses spectral encoding. In 
this method, spreading is achieved by encoding of amplitude 
or phase of the data spectrum [27], [28]. 

This paper is organized as follows: Section [II] provides 
background material along with an overview of the results 



of this paper. In Section III we propose two new bounds 
on the size of 2-D OOCs - we use these bounds to prove 
optimality of some of the constructions presented in the current 
paper as well as of two constructions previously known in the 
literature, but which were not known to be optimal. In the next 
section, we propose five families of 2-D OOCs constructed 
using polynomials over finite fields. In Section [V] we present 
two families of asymptotically optimal codes constructed using 
rational functions over finite fields. We show in Section [Vl] 
how one can generate two asymptotically optimal families by 
composing two of the previous constructions with a constant 



weight binary code. In Section VII we present a three- 
dimensional OOC using polarization as the third dimension. 



In Sections VIII and IX we use generalized bent functions 
to construct two families of efficient asynchronous phase- 
encoding sequences for Optical CDMA. The last section 
concludes the paper. New results are presented as Propositions, 
known results appear as Theorems. Most of the proofs have 
been moved to the Appendices to ensure smooth reading of 
the paper. 

II. Background and Results 



The focus of the entire paper (except for Sections VIII 



and IX i is on direct-sequence encoding. Phase-encoding is 
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TABLE I 

Optimal and asymptotically optimal 2-D optical orthogonal codes in the literature 



Construction 
Name 


Parameters 


Code Size 


Constraint 
Satisfied 


Optimality 


Lee, Seo [29] 


(A x T,3, 1) 


6Ast + As + At, where s and t 
are the sizes of the optimal OOCs 
(A, 3, 1) and (T, 3, 1) respectively. 
{KJ) wnen 1\, 1 =1 ^inoa oj, 
(AO) otherwise 


None 


(O) 


Shurong 
et. al. [30] 


(A x T,£j, 1), 

A — 7)^ 


A(AT-l) 


None 


(O) 


Kwong, 
Yang [31] 


(A x T,lu, 1), 

A = pip2 ■■■pk,T = pi, 

Pk > Pfe-i > • ■ ■ > Pi > ^ 


A "+ a " + a2 +... + A 

Pk PkPk-l PkPk-lPk-2 


AM-OPPTS 


(AO) 


Yang, 


(A x T,u>, 1), 

A — T — dit (tii — "H 4- 1 

A is prime, t is some integer 


A(A 2 -l) 
u;( W -l) 


None 


(O) 


Yang, 

Kwong [32] 


(A x T,lu, 1), 

f/? — A T 1 — T)i -no - - - r> ; 

Pk > Pk-i > ■■■ > Pi > A 


T 


OPPW 


(0)* 


Yang, 

Kwong [32] 


(A x T,w, l),w = A - 1, 

A = pi,T= (pi - l)p 2 ---Pfc, 

Pk > Pfe-i > • ■ ■ > Pi 


AT 
A— 1 


AM-OPPW 


(AO) 


Kwong 
et. al. [33] 


(A x T,lu, 1), 

A = piP2 ■ ■ -Pfc. 

Pfe > Pfc-l > •• • > Pi > u 


A 2 ■ $ooc, 

where ^OOC i s m e cardinality 
of the optimal (T,u, 1) OOC 


AM-OPPTS 


(AO) 


Shivaleela 
et. al. [34] 


(A x T,u>, 1), 

A = T = T is prime 


T 


OPPW 


(0)* 



Here, 

• (O) denotes Optimal, 

• (AO) denotes Asymptotically Optimal, 

• p or pi denotes a prime. 

* These constructions are shown to be optimal using the bounds proposed in this paper. 



restricted to Sections VIII and IX only. 

The advent of Wavelength-Division-Multiplexing (WDM) 
and dense- WDM (D-WDM) technology has made it possible 
to spread in both wavelength and time [32]. The corresponding 
codes are variously called wavelength-time hopping codes and 
multiple- wavelength codes. Here, we will simply refer to these 
codes as two-dimensional OOCs (2-D OOCs). 

A 2-D (A x T,lj, k) OOC C is a family of {0, 1} A x T 
arrays of constant weight u. Every pair {A, B} of arrays in 
C is required to satisfy: 



2-D OOCs 



A T-l 



(1) 



^^A{\t)B{X,{t® T t))<k, 

A=l t=Q 

where either 4 ^ B or r / 0. We will refer to k as the 
maximum collision parameter (MCP) when in addition to ([T| 
holding for all r, we have that equality holds in ([T| for some 
pair A, B and for some r. Note that the asynchronism is 
present only along the time axis. 

Practical considerations often place restrictions on the place- 
ment of pulses within an array. With this in mind, we introduce 



AM-OPPW 

P3, CP1, CR1 



AM-OPPTS 

P4 



No Restriction 



OPPW 

P1, R1 



OPPTS 

P2, P5, R2 



Fig. 1. Various types of 2-D OOCs. The symbols P3, CP1 etc. are reference 
to specific constructions appearing in Table [TT] 



the following terminology (see Fig. [TJ: 

• arrays with one-pulse per wavelength (OPPW): each 
row of every (A x T) code array in C is required to have 
Hamming weight = 1. 

• arrays with at most one-pulse per wavelength (AM- 
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TABLE II 

New 2-D optical orthogonal codes proposed in this paper 



Construction 
Name 


Parameters 


Code Size 


Constraint 
Satisfied 


Optimality 


pi 


(A x T,u,k), 
k < uj = A < T, 
T is prime 




OPPW 


(O) 


P2 


(A x T,u,k), 
A = p, w = T, 

K<CJ,T\p-l 




OPPTS 


(AO) 


P3 


(A x T,u,k), 

1 < A < p m , T = p m — 1, 

CO = A — K, K K UJ 


(T+l)" ; + 1 -l 
T 


AM-OPPW 


(AO) 


P4 


(A x T,ui,k), 

to = T - k, T | p m - 1, 

A = p m - 1, K < UJ 


t Ed|A ((A + i)W 


AM-OPPTS 


(AO) 


P5 


(A x T,ui,k), 
uj = T, A = p m , 

K < UJ, T | p m - 1 


^Ed|(A-l) 


OPPTS 


(AO) 


Rl 


(A x T,u,k), 

ui = A, A < T - 1, 

T = p m + 1, k < is even 


#+1 


OPPW 


(AO) 


R2 


(A x T,uj, «), 

T|p m - 1, A = p m + 1, 

tij = T, k < a; is even 


t( 9 1 -i) Eft(x) u ( d de g(H x )), T , ±)K h ( x )) 

where the sum is over monic h(x) S V q [x\ 
of deg < d 


OPPTS 


(AO) 


CP1 


(A x T,u,k), 

k < a; < A, T is prime 




AM-OPPW 


(AO) 


CR1 


(A x T,u,k), 

uj < A, k < a; is even, 
T = p m + 1 




AM-OPPW 


(AO) 



Here, 

• (O) denotes Optimal and (AO) denotes Asymptotically Optimal, 

• p denotes a prime, q = p m 

• ^t(-) is the Mobius function, 

• /t(-) and u(-) are defined in equations l(6} and j_2| respectively, and 
f q 2t+1 -q, 4 = 1,2,3,4,5,6 

t > 7. 



shift is zero. This is obvious since there is (at most) one 
1 in each row; hence, the time-shifted code matrix cannot 
have any overlap with the original code matrix. We shall use 
this fact later while proving the correlation properties of our 
constructions. 

A. Johnson Bound 

For a given set of values of A, T, ui, n, let $(A x T, u>, k) 
denote the largest possible cardinality of a (A x T,uj,k) 2-D 
OOC. The following adaptation of the Johnson's bound for 
constant weight codes to 2-D OOCs was first noted by Yang 
and Kwong in [32]: 



OPPW): each row of any (A x T) code in C is required 
to have Hamming weight < 1. 

• arrays with one-pulse per time slot (OPPTS): each 
column of every (A x T) code array in C is required 
to have Hamming weight = 1. 

• arrays with at most one-pulse per time slot (AM-OPPTS): 
each column of any (A x T) array in C is required to have 
Hamming weight < 1 . 

The constructions mentioned in Fig. [T] are proposed in this 
paper and have been summarized in Table [TTJ 

Remark 1: Note that for codes that are AM-OPPW or 
OPPW, the autocorrelation for non-zero values of the time 
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Theorem 1: [Johnson Bound] 



$(A x T,u,k) < 



A 



AT - 1 
uj - 1 



AT 



(2) 



A construction meeting this bound with equality is called 
an optimal construction. A construction that meets this bound 
asymptotically (i.e., when A or T (or both) tend to infinity) is 
called an asymptotically optimal construction. 



B. Literature Review 

In this subsection, we focus primarily on prior constructions 
of optimal or asymptotically optimal constructions of 2-D 
OOCs in the literature. A summary of these appears in Table [I] 
Note that all these constructions are optimal (or asymptotically 
optimal) only for MCP = 1. However, the constructions 
that we propose in this paper are optimal (or asymptotically 
optimal) for all values of the MCP, i.e., MCP > 1 thereby 
leading to larger size (this is explained in detail in the next 
subsection). 

The construction by Lee and Seo [29] spreads in the wave- 
length and the time domain by using two different 1-D OOCs. 
Shurong et. al [30] construct a 2-D OOC by employing a fre- 
quency hopping code to spread in the wavelength domain and 
a 1-D OOC to spread along the time axis. The construction by 
Kwong and Yang [31] interchanges the time and wavelength 
components of a frequency-hopping code and then applies 
specific cyclic shifts to control the value of the MCP. The 
first construction by Yang and Kwong [32] uses a 1-D OOC 
to achieve spreading in the wavelength and time domains. The 
remaining two constructions in [32] modify frequency-hopping 
codes to construct 2-D OOCs. The construction by Kwong et. 
al [33] spreads in the wavelength domain using a frequency- 
hopping code and in the time domain using a 1-D OOC. The 
OPPW construction by Shivaleela et. al [34] places a 1 in the 
first time slot of the first wavelength. By cyclically shifting 
the position of the 1 in the subsequent wavelengths by k, the 
entire 2-D code is generated. The k for different codewords 
varies from to T — 1, where T is the number of time slots. 

Additional papers in the literature dealing with the design 
of 2-D OOCs include [35]-[60]. However, since the focus 
of the current paper is on optimal or asymptotically optimal 
constructions, we do not discuss these further here. A paper 
relating to 3-D code construction is [61]. 



C. Overview of Results 

We propose a version of non-binary Johnson bound and 
derive two other bounds from it - these bounds provide upper 
bounds on the size of 2-D OOCs for the case of AM-OPPW 
2-D OOCs and OPPW 2-D OOCs. A special instance is shown 
to lead to the Singleton bound. 

We then propose 9 new families of 2-D OOCs of large 
size. All the codes proposed in this paper are optimal (or 
asymptotically optimal) with respect to the original Johnson 
bound [32] or the new bounds proposed in this paper. We 
obtain codes with large size by constructing optimal families 
for large values of the MCP. Consider MCP = k = 1, for 



example. The 2-D Johnson bound gives 

$(A x T,u,k) 



< 



A 


AT - 1 




UJ 


_ uj- 1 _ 





A /AT - 1 

UJ \ UJ — 1 

A 2 T 



uJ~ 



Similarly, for large values of k, we get 



$(A x T,oj,k) 



< 



K+l 



This shows that, for fixed values of A,T and oj, the 
upper bound on the maximum number of codewords increases 
exponentially in the MCP. Note that all the constructions in 
Table I have MCP = 1 thereby restricting the size to ^pr- 
Hence, by constructing optimal (or asymptotically optimal) 
constructions for larger values of the MCP, we are proposing 2- 
D OOCs with size larger than has been previously constructed. 

All the 9 constructions presented in this paper have been 
summarized in Table \H\ Additionally, Fig. |2]provides a visual 
depiction of the parameter sets for small A, T in the range, 
2<A<17, 2<T<33 for which constructions are now 
available as a result of the 9 constructions presented here. 
All constructions are either optimal or else drawn from a 
family of asymptotically optimal constructions and correspond 
in every case to a code whose size is large in relation to the 
number T of time slots. This table brings out the need for 
proposing different constructions since these 9 constructions 
are applicable to different values of A and T. For ease of 
representation, we use the following legend: 



Construction 
Name 


Name in 
the Table 


Construction 
Name 


Name in 
the Table 


PI 


A 


Rl 


F 


P2 


B 


R2 


G 


P3 


C 


CP1 


H 


P4 


D 


CR1 


I 


P5 


E 







For example, consider the entry in the table corresponding to 
T = 5, A = 3. The entry reads AFHI. This means that for T = 
5 and A = 3, we are proposing four optimal (or asymptotically 
optimal) constructions in this paper, viz., construction PI, Rl, 
CP1 and CR1. 

Remark 2: Note that the rows corresponding to T = 21, 25 
and 27 are empty. However, this does imply that there are 
no constructions for these value of T - it simply means that 
there are no constructions for 2 < A < 17. For example, for 
T = 21, we have construction P2 ( or B) for A = 43, 127, 211 
etc. Similarly, we have construction P2 for T = 25, A = 
101,151,251; and for T= 27, A = 109,163,271 and so on. 

The seven constructions (PI to P5, Rl and R2) are generated 
by regarding a codeword in the 2-D code as the graph of 
a function. For example, a codeword in a 2-D AM-OPPW 
code can be regarded as the graph of a function t = /(A), 
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< t < T — 1, < A < A — 1 mapping wavelength to time. 
Analogously, a codeword in a 2-D AM-OPPTS OOC can be 
regarded as the graph of a function X = f(i) , < t < T — 1, 
< A < A — 1 mapping time to wavelength. 

The first 5 of these codes (PI to P5) are constructed using 
polynomials over finite fields. Code PI is optimal while codes 
P2 to P5 are asymptotically optimal. 

The next two codes (Rl and R2) are constructed using 
rational functions over finite fields. Both these codes are 
defined for even values of the MCP and are asymptotically 
optimal. 

Next, we compose PI with a constant weight binary code 
and generate the concatenated code CP1. This code is asymp- 
totically optimal and is AM-OPPW. We do a similar compo- 
sition of Rl and a constant weight binary code to generate 
CR1, which is also asymptotically optimal. 

We present a 3-D code construction using polarization as the 
third dimension. This code is constructed using the Chinese 
Remainder Theorem. 

Finally, we use generalized bent functions to construct a 
family of efficient asynchronous phase-encoding sequences for 
optical CDMA. 

III. New Bounds on the Code Size 

In this section, we begin by proposing a one-dimensional 
Johnson bound on constant weight codes over a non-binary 
alphabet. The bound will establish the optimality of some of 
the constructions that we propose in the sections that follow. 

The codes under consideration are over an alphabet A of 
size (T + 1) containing a distinguished element, which we 
shall call 0. For codes over such an alphabet, we define 
the Hamming correlation between two codewords to be the 
number of symbol locations in which the two codewords 
contain the same non-zero symbol. 

Let At(A,ui, k) denote the maximum possible size of a 
constant-weight code C over the alphabet A of size T + 1 of 
length A, Hamming weight u>, and Hamming correlation < k. 



Proposition 2 (Nonbinary Johnson Bound): 



A t (A,cj, k) < 



TA 

to 



T(A - 1) 
uj - 1 



T(A - k) 



Proposition 3 (Bound on AM-OPPW Code Size): The 
maximum possible size <1>(A x T, u>, n) of an AM-OPPW 2-D 
OOC with parameters (A x T,u>, k) satisfies: 



$(A x T,lu,k) < 



T(A - 1) 



T(A - k) 



For the proof, we refer the reader to Appendix [I] 



For the proof, we refer the reader to Appendix [TTJ 

Corollary 4 (Bound on OPPW Code Size): The maximum 
possible size $(A x T,u,k) of an OPPW 2-D OOC with 
parameters (A x T,lo,k) satisfies: 

$(A x T,u,k) < T K . 
It follows from the above that the following constructions 
are optimal (note that k = 1 here) although this was unknown 
to the authors (see Table [[}: 

• the second construction by Yang and Kwong [32] and 

• the construction by Shivaleela et. al [34]. 

The next few sections will each introduce new constructions 
of optimal (or asymptotically optimal) 2-D OOCs. 

IV. New Constructions Based on Polynomial 
Functions 

A codeword in a 2-D AM-OPPW OOC can be regarded as 
the graph of a function f = /(A), < t < T-l, < A < A-l 
mapping wavelength to time. Analogously, a codeword in a 2- 
D AM-OPPTS OOC can be regarded as the graph of a function 
A = /(f), < f < T - 1, < A < A - 1 mapping time to 
wavelength. 

Without loss of generality, in the constructions below, we 
will identify the T time slots with the set Zy, which we label 
as T. This will allow us to identify cyclic shifts in the time 
domain with (mod T) additions in Z^ . We will identify the 
A wavelengths with subsets C of algebraic structures. Thus, 
| T |= T and | £ |= A. 

All of the constructions in this section employ polynomial 
functions whose degree is bounded above by the desired value 
of the MCP k. In the next section, we use rational functions 
to construct 2-D OOCs. 

While the constructions below have some elements in com- 
mon, they also have their differences. The reason for providing 
the set of constructions is to provide constructions for as many 
(A, T) parameter sets as possible. The impact of the different 
constructions given here can be seen in Fig. [2j where we 
identify constructions with the entries in the table. 



Remark 3 (Recovering the Binary Johnson Bound): In the 
binary case, i.e., when T+l = 2, the above inequality reduces 
to the Johnson bound for constant-weight binary codes. 

Remark 4 (Singleton Bound): For the special case when 
A = uj, the above bound reduces to 

A t {A,u,k) < T K+1 

and we have, in fact, recovered the Singleton bound of coding 
theory. 

We now proceed to apply the non-binary Johnson bound to 
derive bounds on AM-OPPW 2-D OOCs. 



A. Construction PI: Mapping Wavelength to Time, OPPW, 
u) = A, A < T, K < uj, T prime 

Here T = p, T = Z p , C C Z p and we consider 
polynomials /(A) over Z p of degree < n mapping C — * T, 
where p is a prime. Let (p be a mapping from the set of all 
such functions to the code matrices, where the A x T code 
matrix C associated to function / is given by C(A, t ) = 1 iff 
/(A) = t, where A e C. 

We first show that tp is an injective mapping, i.e., the code 
matrices C\ and C2, corresponding to polynomials fx and 
/2 respectively, are equal iff /1 = fa. This is true since 
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CEFI 


ACH 


CFI 


iZ 


CFI 


AH 


CFI 


AH 


LL 


o 


o 


AFHI 


CFI 


AH 


\L 




o 


AH 


CFI 




CFI 




CFI 


AH 


CFI 


AH 


iZ 


\L 






DGH 


CHI 


CDGI 


AEFHI 


CFI 


ACH 


CFI 


LL 


CFI 


AH 


CFI 


AH 


LL 


o 


u 


AFHI 


CFI 


AH 


LL 




o 


AH 


CFI 




CFI 




CFI 


AH 


CFI 


AH 


iZ 


\L 




CO 


BCEH 


ACDGHI 


CEFI 


AFHI 


CFI 


ACH 


CFI 


LL 


CFI 


AH 


CFI 


AH 


LL 


o 


o 


AFHI 


CFI 


AH 


iZ 




o 


AH 


CFI 




CFI 




CFI 


AH 


CFI 


AH 


iZ 


UL 




CN 


ACDH 


ACEFHI 


CFI 


AFHI 


CFI 


ACH 


CFI 


iZ 


CFI 


AH 


CFI 


AH 


LL 


o 


o 


AFHI 


CFI 


AH 


iZ 




o 


AH 


CFI 




CFI 




CFI 


< 


CFI 


AH 


iZ 


UL 




T \ A 


CM 


CO 




LO 


CD 




00 




o 




CM 


CO 




LO 


CD 




00 




O 
CM 


CM 


CM 
CM 


CO 
CM 


^t- 

CM 


LO 
CM 


CD 
CM 


h~ 
CM 


00 
CM 


CM 


o 

CO 


CO 


CM 
CO 


CO 
CO 






I- 



Fig. 2. Representative table showing constructions for (A, T), A G [2, 17], T S [2, 33] 
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d(X,t) = C 2 (X,t) A(A) = / 2 (A). Since /i(A)-/ 2 (A) is 
a non-zero polynomial of degree < k, this equation can have a 
maximum of k zeroes. Hence, as ui > k , these two functions 
cannot coincide in all the to positions. 

Since this code is OPPW, the autocorrelation function for 
each of these matrices for non-zero cyclic shifts is by 
Remark Q] 

We next declare two polynomials /(A) and /(A) +5, where 
8 E Z p , to be equivalent. This causes the set of all the code 
matrices to be partitioned into equivalence classes, each of 
size T. This results in T = T K code matrices. 

Note first that C a (X,t + r) = 1 = C b (X,t) iff / (A) - 
t = /b(r). For cross-correlation, consider C a (X,t + r) and 
Cfc(A,i) from different equivalence classes corresponding to 
functions f a and respectively. Next, the polynomial f a (X) — 
t — fb{X) is non-zero since f a and fb belong to different 
equivalence classes. Since the polynomial is non-zero, it can 
have a maximum of k zeroes (the degree of the polynomial); 
hence, the two code matrices C a and Cf, can have a maximum 
of k collisions. 

This results in a (A x T, A, k) 2-D OOC of size T K . This 
construction can be seen to be optimal by Corollary [4] 

Remark 5: The elements of this 2D-OOC can also be 
regarded as corresponding to codewords in a Reed-Solomon 
code under an appropriate equivalence relation between the 
codewords. Let A,k,T be as above. In particular T = p, 
where p is a prime. The [A, k] Reed-Solomon code Crs ma Y 
be constructed as follows: let {ai, a 2 , ■ ■ ■ , a\} denote a set 
of A distinct elements drawn from ¥ p . Let V K denote the set 
of all polynomials over F p of degree < n. Set 

Crs = {(f(ax),...,f(a A )) \feV K }. (3) 

Next, partition the set of all the codewords into p K equivalence 
classes by declaring c x ~ c 2 if c 1 — c 2 = i]l , r\ S F*. Finally, 
we form a set S by picking precisely one element from each 
equivalence class and by associating a (A x T) matrix A(X, t) 
to each vector a E S by setting 

A(X, t) = 1 iff a x = t, 1 < A < A. 

B. Construction P2: Mapping Time to Wavelength, OPPTS, 
A = p, p prime, T \ p — 1, u = T, k < to 

Here C = Z p . Let a be an element of Z p of multiplicative 
order T and let H be the subgroup of order T generated by a 
in Z p . We will identify 7Lt with H by associating t with a'. 
Consider polynomials f(x) over Z p of degree < k mapping 
H — * C. Let (p be a mapping from the set of all such functions 
to the code matrices, where the A x T code array C is obtained 
by setting C(X,t) = 1 iff /(a 4 ) = A, where t <E T and 
A G Z p . 

We first prove that is injective, i.e., the code matrices C\ 
and C-2, corresponding to polynomials fx and f 2 respectively, 
are equal iff /j = f 2 . This is true since Ci(A, t) = C 2 (X, t) O 
/i(a*) = f 2 (a f ). Since fi(a f ) — f 2 {a l ) is a non-zero 
polynomial of degree < k, this equation can have a maximum 
of k zeroes. Hence, as u> > n, these two polynomials cannot 
coincide in all the u> positions. 



We next discard all sub-period polynomials, i.e., polynomi- 
als f{x) that satisfy /(a z x) = fix) for some i 6 V T . This 
ensures good autocorrelation. Since fia l x) — fix) is not the 
zero polynomial, we know that it has a maximum of k zeroes 
(the degree of the polynomial). Hence, the autocorrelation is 
bounded above by k. 

We now define two polynomials fix),gix) to be equivalent 
if fiofx) = gix) for some i E Zy- We pick a code matrix 
corresponding to each equivalence class to form a code of size 
T X/d|(A-i) LAJt-I — ij nid), where /i(-) is the Mobius 
function (see [5] for code size computation). 

Now, consider two polynomials f a and drawn from 
different equivalence classes. Consider the corresponding code 
matrices C a and Cfc. We know that C a (A, t + r) — Cf,(A, t) 
/a(a t+T ) = /&(«*)■ Since these two polynomials have been 
drawn from distinct equivalence classes, the difference polyno- 
mial faict T x) — fbix) is non-zero. This implies that f a i& T x) — 
fbi%) has < k zeroes. This proves that the crosscorrelation is 
bounded above by k. 

This construction is shown to be asymptotically optimal in 
Appendix [Hi] 

C. Construction P3: Mapping Wavelength to Time, AM- 
OPPW, u = A- k, K <lu, 1 < A < p m , T = p m - 1, p 
prime 

Let p be prime and a a primitive element of F p m. Let 
T = p' m — 1, T = Zt and L C F p n,. Consider polynomials 
fix) over Fpm of degree < k mapping C —> F* m . Let ip 
be a mapping from the set of all such functions to the code 
matrices, where the A x T code array C is obtained by setting 
C(A, t) = 1 iff fiX) = a\ where t E T and A G C. For those 
values of A such that /(A) = 0, the entire row is left blank. 
Clearly, at most k rows in any matrix can be blank. To restore 
the constant weight property, we arbitrarily delete appropriate 
number of l's to keep the weight equal to A — k for all the 
codewords. 

We first prove that cp is injective, i.e., the code matrices C\ 
and C 2 , corresponding to polynomials fi and f 2 respectively, 
are equal iff /j = f 2 . This is true since Ci(A, t) — C 2 iX, t) O 
/i(A) — f 2 iX). Since /i(A) — fiiX) is a non-zero polynomial 
of degree < k, this equation can have a maximum of k zeroes. 
Hence, as lu > k, these two polynomials cannot coincide in 
all the lu positions. 

Since this code is OPPW, the autocorrelation function for 
each of these matrices for non-zero cyclic shifts is by 
Remark Q] 

We define two polynomials fix), gix) of degree < k 
to be equivalent if a 1 fix) — gix) for some i E Zt- 
In our construction, we choose one polynomial from each 
equivalence class. This gives us a code of size ( - T+1 i r — — . 

Now, consider two polynomials f a and fb drawn from 
different equivalence classes. Consider the corresponding code 
matrices C a and Cb- We know that C a (A, t + r) = Cb(A, t) <=> 
a ~ T fai x ) — fbi%)- Since these two polynomials have been 
drawn from distinct equivalence classes, the difference poly- 
nomial oT T faix) — fbix) is non-zero. This implies that 
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a ~ T fa{x) — fb{x) has < k zeroes. This proves that the 
crosscorrelation is bounded above by k. 

This construction is proved to be asymptotically optimal in 



Appendix IV 



D. Construction P4: Mapping Time to Wavelength, AM- 
OPPTS, A = p m -l, p prime, uj = T-n, k < w, T | p m -l 

Let p be prime. Let a be a primitive element of F pm and 
let (3 be a non-zero element in ¥ pm of multiplicative order 
T. Let H be the subgroup of F p ™ generated by (3. We will 
identify Z T with H by associating t with /3*. Let C = F* m = 
{1, a, . . . ,a p _2 } and T = Zt. Consider polynomials f(x) 
over F pm of degree < n mapping H C. Let ip be a mapping 
from the set of all such functions to the code matrices, where 
the A x T code array C is obtained by setting C(X,t) = 1 
iff /(/3*) = A. For those values of A such that /(/?') = 0, the 
entire column is left blank. Clearly at most k columns in any 
matrix can be blank. To restore the constant weight property 
we arbitrarily delete an appropriate number of l's to keep the 
weight equal to T — n for all the codewords. 

We first prove that p is injective, i.e., the code matrices C\ 
and C 2 , corresponding to polynomials fx and / 2 respectively, 
are equal iff fx = ,f 2 . This is true since Ci(A, t) — C 2 (A, t) 

= / 2 (/3*). Since A(/3*) - / 2 (/?') is a non-zero poly- 
nomial of degree < k, this equation can have a maximum of 
k zeroes. Hence, as ui > k , these two polynomials cannot 
coincide in all the u> positions. 

We next discard all sub-period polynomials, i.e., polynomi- 
als f(x) that satisfy f(f3 l x) = f(x) for some i G Z^,. This 
ensures good autocorrelation. Since f(f3 l x) — f(x) is not the 
zero polynomial, we know that it has a maximum of k zeroes 
(the degree of the polynomial). Hence, the autocorrelation is 
bounded above by k. 

We now define two polynomials f(x),g(x) to be equivalent 
if f(f3 l x) = g(x) for some i G Zt- We pick a code matrix 
corresponding to each equivalence class to form a code of 
size J2d\A f (A + 1) — ij fi(d) (see [5] for code size 
computation). 

Now, consider two polynomials f a and fb drawn from 
different equivalence classes. Consider the corresponding code 
matrices C a and Cb- We know that C a (X,t + r) = Cb(X,t) 
fa(f3 t+T ) — fb (/?*)■ Since these two polynomials have been 
drawn from distinct equivalence classes, the difference polyno- 
mial f a (p T x) — fb(x) is non-zero. This implies that f a (f3 T x) — 
fb(x) has < k zeroes. This proves that the crosscorrelation is 
bounded above by k. 

Since the number of codewords is similar to the number 
of codewords for construction P2, the proof for asymptotic 
optimality is along the same lines - see Appendix III 



E. Construction P5: Mapping Time to Wavelength, OPPTS, 
A = p m , p prime, lu = T, k < uj, T\ p m - 1 

Let p be prime. Let a be a primitive element of ¥ pm 
and let (3 be a non-zero element in F p ™ of multiplicative 
order T. Let H be the subgroup of F p ™ generated by (3. 
We will identify Zt with H by associating t with (3*. Let 



C = F p m = {0,1, a, ...,a p ~ 2 } and T = Z T . Consider 
polynomials f(x) over F p m of degree < k mapping H — * C. 
Let ip be a mapping from the set of all such functions to the 
code matrices, where the A x T code array C is obtained by 
setting C(X,t) — 1 iff f((3 f ) — A. Thus, the construction here 
is along the lines of the previous construction except that the 
wavelengths are in 1-1 correspondence with all of ¥ pm . 

We first prove that tp is injective, i.e., the code matrices C\ 
and Ci, corresponding to polynomials fx and / 2 respectively, 
are equal iff fx — / 2 - This is true since Ci(A, t) = C 2 (A, i) 
/i(/3 4 ) = / 2 (/3'). Since - / 2 (/3 4 ) is a non-zero poly- 

nomial of degree < k, this equation can have a maximum of 
k zeroes. Hence, as u> > k , these two polynomials cannot 
coincide in all the uj positions. 

We next discard all sub-period polynomials, i.e., polynomi- 
als f(x) that satisfy f(f3 l x) = f(x) for some i G Z* T . This 
ensures good autocorrelation. Since f(f3 % x) — f(x) is not the 
zero polynomial, we know that it has a maximum of k zeroes 
(the degree of the polynomial). Hence, the autocorrelation is 
bounded above by k. 

We now define two polynomials f(x),g(x) to be equivalent 
if f(f3 l x) = g(x) for some i G Zt- We pick a code matrix 
corresponding to each equivalence class to form a code of 

size ?£d|(A-l) ( A ^l 
computation). 

Now, consider two polynomials f a and drawn from 
different equivalence classes. Consider the corresponding code 
matrices C a and Cb- We know that C a (X, t + r) = Cb(X, t) 
fa(P t+T ) = /f)(/?*). Since these two polynomials have been 
drawn from distinct equivalence classes, the difference polyno- 
mial f a ((3 T x) — fb(x) is non-zero. This implies that f a (f3 T x) — 
fb(x) has < k zeroes. This proves that the crosscorrelation is 
bounded above by k. 

The proof for asymptotic optimality of construction P5 is 
similar to the proof for construction P2 - see Appendix III 



1 ) (see [5] for code size 



V. New Asymptotically Optimal Constructions 
Based on Rational Functions 

We now use rational functions over F g to construct two 
2-D OOCs. Some properties of rational functions that we 
use are summarized in the first subsection; the following two 
subsections deal with the two constructions. 

A. Cyclic Ordering for the Projective Line 

Consider all elements of F 2 other than the element [0 0] T . 
We define an equivalence relation amongst these elements as 
follows: 



(4) 

-qb] T . This 



[a b] T - [c df, 

if for some ?/ G F*, we have [c d] T = [r/a 
partitions F 2 into (q+ 1) equivalence classes, with each class 
containing (q — 1) elements. The projective line is obtained 
by taking precisely one element from each equivalence class. 
We denote this by P 1 (F g ). Thus, there are (q+ 1) "points" on 
the projective line. We will use 

a 
b 
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to denote the equivalence class containing 



We next present a cyclic ordering of the elements of P 1 (F g ). 
Let h(x) = x 2 + h\X + h with hi, ho € ¥ q be a primitive 
polynomial over ¥ q . Let a be a zero of this polynomial and 



H 



-ho 

1 -hi 



be the associated companion matrix. 
Theorem 5: We claim that 



H 1 



eq 



< i < q 



(5) 



Proof: Note from the definition of the companion matrix 



that 



r 




a 


o_ 




P. 



is equivalent to saying that 

a 1 ■ 1 = a + ba. 



Thus 



T 




T 




0_ 








with j > i iff 

a 1 = 6a j 

for some 6> e F*, i.e., iff 

(9 + 1) I U-i)- 
Hence, the equivalence classes 



H' 



eq 



< i < q 



are all distinct, thus proving the theorem. 

■ 

Example 1: We present a cyclic ordering of P 1 (F 3 ). The 
polynomial f(x) = x 2 + x + 2 is a primitive polynomial over 
F 3 . Thus 



H 
H 





" 




-2 " 






" 


1 " 










1 




-1 






1 


2 








1 











1 




1 




1 




2 




. H" 















1 











2 









2 



1 









1 




2 





eq 


i — 1 


eq 


_2_ 


eq 


_2_ 



P^Fs) = 
It can easily be checked that 

which shows that the ordering is cyclic. 



eq, 



1 




2 




1 — 1 
















B. Construction Rl: Mapping Wavelength to Time, OPPW, 
T = q + 1, for q = p m , p prime, A < q, uj = A, n < v, 

k = 2k! = 2d 



Set T = P^F,) and L 



Let /, g be polynomials over ¥ q and let 
mapping given by: 



{ai,a 2 , . ■ . ,a\, . . . ,(Xa} C 
denote the 



,9(«i) 



1 < i < A. 



eq 



We will regard as a rational function map because when 
g(cti) ^ for any i, we equivalently have: 



(j)(ai) 



Ha) 
g(<*i) 
1 



eq 



Let Td denote the class of rational functions 

~m 

over F q (and hence over C by restriction) where / and g are 
polynomials over F g satisfying: 

1) / and g are both non-zero and both of degree < d, 

2) f and g are relatively prime, i.e., (/, g) = 1, 

3) / is monic, and 

4) either / or g must be a non-constant function; equiva- 
lently, ^ifv is not the constant function. 

9\ x ) 

The last condition has been included here since it is mathe- 
matically convenient to exclude the constant functions at this 
stage and bring them back later. Let //(•), ji(-) : F q [x] — ► Z 
be the function defined by 



£(&(*)) 



n, &(s) = i. 

(— l) r , b(x) is the product of r monic, 

irreducible polynomials over ¥ q , 
0, else. 



(6) 



Then the number Cd =\ Td \ of rational functions in Td can 
be computed from the results in [5] and is given by: 

(qds+i _ 1)2 

Jq^W 



T d 



E 

h(x) 



-KKx)) - (q-l) 



where the sum is over all monic polynomials h(x) e F q [x] 
of degree s < d and where the last term accounts for the 
disallowed constant functions. It can be shown [5] that 



-,2d+l 



C(d) 



q, 



> g 2d+i _ r 



d= 1,2,3,4,5,6 
d > 7. 



(7) 



Let <p be a mapping from Td to the code matrices, where 
the A x T code array is obtained by setting C(A, t) = 1 iff 



0("a). 



We first prove that <p is injective, i.e., the code matrices 
C a and Cb, corresponding to rational functions and 

[ll(a X x )\ respectively, are equal iff f a (x) = f b (x) and g a (x) = 
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gb(x). Clearly, the code matrices C a and Cb are equal iff for 
every 1 < A < A, 



'fa(ax)' 




'fb(ax)' 









for some 9\ 6 F* 



(8) 



i.e., iff 



6\gb{ot\) V A 
OVA 
OVA 

as polynomials 

fb(x)g a (x). 



fa{a\) = 9\fb(a\) and g a (a x ) 
& 8\ [fa(a\)gb(a x ) - fb(u\)g a (a\)} 
^ fa{u\)gb{u\) - fb(a\)g a (a\) 
fa{x)gb{x) - f b (x)g a (x) 
<=> fa(x)g b (x) 

Since / a (x) is coprime to g a (x), it must be that f a (x)\fb(x). 
Similarly, f b (x)\f a (x). Since both f a (x) and fb(x) are monic, 
it implies that f a (x) = fb(x) which forces g a (x) = gb(x) and 
hence proves that (p is injective. 

Since this code is OPPW, the autocorrelation function for 
each of these matrices for non-zero cyclic shifts is by 
Remark Q] 

Next, we would like to establish that the collection of 
matrices ip(Td) is closed under cyclic shifts along the time 
axis. To show this, it is sufficient to show that if C(X,t) is 
the code matrix associated to element [^j] € Td, then 



there exists 



17 (»)] 

Lg'(x)J 



C(X,t+l (mod T)) 
found by setting 



Td whose associated code matrix is 
Such an element | ' , )*{ | 6 Td is easily 



[/ (x)-l 



'f'(x)- 




' -ho ' 




'fix)' 




m'( x ). 


= Mo 


1 -ht 






= Mo 



-h Q g(x) 
f(x) - hxg(x) 

where /iq is chosen to ensure that (mo)( — ho)g(x) is monic. 

Let us define two elements T^Hl and P^Hl to be equiv 
alent, i.e., 



'fa(xy 




'fb(x)' 


_9a(x)_ 




.9b( x ). 



if the corresponding code matrices C a and C& obtained via 
the mapping ip, are cyclic shifts of each other. Since ip is a 
1 — 1 mapping and since every code matrix C(A, t) in the 
image of <p has T distinct cyclic shifts, it follows that there 
are T elements within each such equivalence class. Our 2D- 
OCDMA code C is then the code obtained by selecting one 
element of Td from each equivalence class and then applying 
the mapping <p. It follows then that this code is of size S^p+l, 
where the 1 arises from the inclusion of the constant function 
and where c(-) is as defined in |7]). 

It remains to establish the cross-correlation bound for the 
code C. Consider two rational functions [^"f^l, \^ h ) x I] and 

their corresponding code matrices C a and Cb- Let [g^]] ^ e 
the element of Td associated to a cyclic shift of C a (X,t) by 
t. Let the value of the cross-correlation be v and let V denote 
the collection of (t, A) € {1, 2, • • • , A} x {1, 2, • • • , T} such 
that 

C a (X,t + T (modT)) = C b (X,t) = 1. 
This implies that for every A such that (A, t) € V, 



fc(a\)' 
.9c( a \). 



fb{a\) 
.9b{a\). 



for some 6\ e ¥* 



(9) 



i.e., 

/cK) = 0\fb(a\) and g c {a\) = 6\g b {a\) 

& Q\ [fc{ot\)gb{u\) - fb(a\)g c (a\)] = 0, 

fc(a\)gb(a\) - fb{ax)g c {a\) = 0. 

But since the product polynomials f c (x)gb(x) and fb(x)g c (x) 
both have degree that is bounded above by 2d = k, it follows 
that | V |< k, thereby establishing that all cross-correlations 
for any cyclic shift are bounded above by k, i.e., that the 
MCP < k and we are done. This construction is proved to be 
asymptotically optimal in Appendix |V] 

C. Construction R2: Mapping Time to Wavelength, OPPTS, 
lu = T, T | q — 1, A = q + 1, for q = p m , p prime, 
k = 2k' = 2d < uj 

We associate with each wavelength a distinct element of 
P 1 (F (J ). Unlike in the case of Construction Rl, the ordering 
of elements of P 1 (F g ) is immaterial here. We associate /3* 
with the t-th time slot, < t < T - 1, where f3 e ¥ q has 
multiplicative order T, T\(q — 1). As before, we associate 
elements along the time axis with elements of Zt- 

We restrict our rational functions to once again belong to 
the set Td for d — | , where Td is as defined in the previous 
subsection. In this case, we need to necessarily exclude the 
constant functions since they lead to constant autocorrelation 
function equal to uj for all time shifts and uj > k. Thus unlike 
the case of Construction Rl, the last condition in the definition 
of Td is needed here. 

Let <p be a mapping from Td to the code matrices, where 
the A x T code array associated to [gSj] is obtained by setting 
C(A,t) = liff 



A. 









3a (P). 


= 0t 





where t e Z T and A e P x (F g ). 

Once again, the first objective is to prove that cp is injective, 
i.e., that the code matrices C a and C&, corresponding to 
rational functions and [^[5] respectively, are equal 

iff f a (x) = fb(x) and g a {x) = g b (x). Clearly, C a and C b are 
equal iff 



where 9 t € F! all < t < T - 1. 



This equation is similar to equation ([8]) obtained in the 
previous subsection and hence by arguing in similar fashion, 
we conclude that ip is injective. 

If C(X, t) is the code matrix corresponding to rational 
function [^j] , and p,f is the coefficient of the highest degree 
term in }{j3 T x), then it is clear that C(X,t + r (mod T)) 
is uniquely associated to — P!l x )] . Next, note that if the 
autocorrelation function of C(A, t) is equal to uj for some 
nonzero time shift r, this must mean that the matrices C(A, t) 
and C(A, t + r (mod T)) are identical, which in turn can only 
happen by the 1 — 1 nature of ip, if and only if the functions 
and —\ f( inl x }] are identical. It follows then that to 
avoid constructing code matrices C(X,t) with autocorrelation 



11 



function equal to lo for some nonzero time shift r, we need 
to discard sub-period functions, i.e., functions that satisfy 





= e 




_g(P T x)_ 







e f* 



(10) 



for any 1 < r < T — 1. Accordingly, our next objective is 
to calculate the size of the set resulting from discarding the 
sub-period functions belonging to Let us first calculate the 
number of polynomial pairs (f(x),g(x)) satisfying 

1) / and g are both non-zero and both of degree < d, 

2) f{x) is monic and 

3) the pair (f(x) 7 g(x)) is not sub-periodic, i.e., that [^S] 
does not satisfy ( |T0] > for any value of time-shift pa- 
rameter r (this condition subsumes the requirement of 
elements of Td, that either f(x) or g(x) be a non- 
constant polynomial). 

We will subsequently modify this count to ensure that [g^j] 
satisfies the remaining requirement that will cause [ to 



belong to namely that 

4) (f,g) = l. 

Let 

r 

f(x)=Y,fiX ei , 

i=l 

where without loss of generality < e% < e 2 < 
and f r = 1. Similarly, let 



< e r < d 



3=1 



9j- 



where < a\ < a 2 < . . . < a s < d. It is straightforward to 
show that (fTOb holds iff for some 1 < r < (T - 1): 



^r(ei-ei) _ 
^r(ei-ai) _ 



1 for 2 < i < r, 
1 for 2 < j < s, 
1. 



and 



(11) 



It will be found convenient to partition the different ways in 
which this can happen as follows: 

• Case (A) e\ = a\ and r = s = 1. 

• Case (B) e\ = ay, either r > 1 or s > 1 and 

gcd ({e, - ei}[ =2 , {a, - ai}^ =2 , T) 

= I for some I > 1. (12) 

. Case (C) ei > ai and 

gcd ({e, - ei}- =2 > - a i}']=2, ( e i - T ) 

= I for some I > 1. (13) 

. Case (D) ai > ei and 

gcd ({e, - ei}- =2 , {aj - ai}* =2 , (ai - ei), T) 

= ^ for some I > 1. (14) 

Our interest is in counting the number of rational functions 
[g(x)] corres P onc ling to Cases B, C, D where I = 1 in each 
case and will carry out the count for the different cases 
separately. Clearly the count for Case C and Case D is the 



same, so it suffices to obtain a count for Cases B and C. We 
begin with the count for Case C. 

Case C Count For l\T, let us define uc(l) to be the number 
of integer sets corresponding to Case C, i.e., the number of 
integer sets, {ei \ < < d, 1 < i < r}, {a,j \ < ay < 
d, 1 < j < s} where ei > a x and where in addition, 



gcd ({e, - ei}^ =2 , {aj - ai}^ =2 , |d - ai|, T) 



I. 



(16) 



Our interest lies in computing wc(l) an d we will do this using 
Mobius inversion [5]. 

Let yc(l) be the number of integer sets {e^ | < ei < 
d,l < i < r}, {aj \ < aj < d, 1 < j < s} where e\ > a\ 
and where in addition, 

. l\T, 

. l\(ei ~ ei) V 2 < i < r, 
• l\(cij — ay) V 2 < j < s and 
. l\(ei - ay) . 
It follows that 



y c (l)= E u c {l'). 

l': l\l'\T 



(17) 



Set 



yc{l) = yc(j), 



T 

■V 



uciP) = u c(j)- 



Then we can rewrite ( fT7| ) in the form 



Z : i|;'|T 



i':i'|T and 



T 



i.e., for any divisor I of T, 



£ C (o = E^(z')- 

Our goal is to compute uc(l) = uc(T). Using Mobius 
inversion allows us to write 



m c (1) = u c {T) 



where ^i(-) is the Mobius function. The value of j/c(0 can be 
shown to be given by 

Vc{l) = 

{^ J+1 -l}E{^ Lft£f±£lJ+1 -l}- d9) 



d 



ei=0 

Since we are interested in counting the number of polynomials 
[g(x)] wne re f(x) is monic, we are interested in the quantity 
J°Sxl which can be computed by substituting for yc from ( |T9] > 
into (ff8|l and then dividing by (q — 1). 
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\c\ 



1 y 

d-—deg(h{x)) 



d—deg(h(x)) 



£ 2 S I 

t|T I ei=0 



|^ d-deg(h(a:))-e 1 j 



|^ d-deg(h.(a;))-e 1 +ct j 



e!=0 



£ ^--i} -(,-i)'(j-*#(>))+i) U(<) 



jl{h{x)). 



(15) 



Case B Count An analogous argument can be used to show 
that the corresponding expressions for ysiP), in Case 

B are given by 



Vb{1) = 

Mi) = 



d 

E 

ei=0 



{^ J+1 -l} 2 -(?-l) 2 (d+l), 



$>b(0a*(0 (20) 

in which the subtracted second term in the expression for 
Ub(1) ensures that the instances when both f(x) and g(x) 
are monomials of the same degree are not counted. 

Overall Count Putting together Cases B, C, D, we see that 
the analogous expressions for y(l) and u(l) for the desired 
overall count u(l) are given by 



y{i) 



ei=0 



2 ^{ (Z L^j + i_ 1 }^{ (Z L^j + i_ 1 } 

c=l 

XI {g L ^ J+1 - i} 2 - - i) 2 + i) , 



(21) 



ei=0 



U(l) 



E 

l\T 



y(l)fi(l). (22) 



The overall count is computed by substituting for y from plj ) 
into ( p2| and then dividing by (q — 1). 

We next proceed to modify this count to ensure that f(x) 
and g(x) are co-prime. 

Let us set 

where we have written u(d, T, 1) in place of w(l) to emphasize 
that u is a function of both d, T as well. In the ensuing, it will 
be found convenient to keep in mind that for any polynomial 



h{x), the function 



rf(x)h(x)i 

lg{x)h(x)\ 



is sub-periodic iff the function 



[g(xj] ' s su b-P er i°dic. Let h(x) be a fixed, monic polynomial 
over ¥ q of degree s and let M(d, T, h(x)) denote the number 
of polynomials satisfying 

1) / and g are both non-zero and both of degree < d, 

2) f(x) is monic 

3) f{x) = h(x)f'(x), g(x) = h(x)g'(x), for some f'(x), 

4) (/ {x),g (x)) are not sub-periodic, i.e., that [yS] does 
not satisfy ( fTO] ) for any value of time shift parameter r. 

It is not hard to show that M(d, T, h(x)) is a function of the 
polynomial h(x) only through its degree s and that moreover, 

M{d,T,h(x)) = N(d-s,T). 



Let us also define M(d, T) to be the set of all polynomials 
[Jjxj] satisfying 

1) / and g are both non-zero and both of degree < d, 

2) f(x) is monic 

3) (f(x),g(x)) are not sub-periodic, i.e., that [^^|] does 
not satisfy ( fTO] ) for any value of time shift parameter r 
and in addition 

4) (/,<?) = ! ■ 

Then we can show that M(cZ, T) is given by 

M(d,T) = ^M^T.AW)^)), 

h(x) 

= J2N(d-deg(h(x)),T)fi(h(x)), 

h{x) 

where the function jl is as defined in |6| and where the 
sum is over all monic polynomials h(x) E F 9 [x] of degree 
< d. We can now define equivalence classes on the M(d, T) 
functions that remain as before by defining the functions 

[ff(xj] anc * /I7 [g(|3 T K)] to ^ e ec l u i va l en,; an d choosing precisely 
one function from each equivalence class. Since there are no 
sub-period functions, there are precisely T elements in each 
equivalence class and thus the total number of code matrices 
is then finally given by 



M[d,T) 



(23) 



and we have completed our count. This overall count is given 
in {[5]). 

It is relatively easy to verify that the autocorrelation function 
is bounded above by 2d = k for all nonzero time shifts 
and that the crosscorrelation function is uniformly bounded 
above by n as well. The proof for asymptotic optimality for 
construction R2 is given in Appendix |VI] 

VI. New Asymptotically Optimal Constructions 
based on Concatenation 

Consider a (AxT, ui, k) 2-D OOC C under the requirement 
that there is at most one pulse per wavelength (AM-OPPW). 
We present two asymptotically optimal constructions that are, 
in a sense, concatenation of a constant-weight binary code 
and an OPPW code. We use this method to construct two new 
AM-OPPW codes by composing a constant weight binary code 
with code PI (or Rl). 

Let C cw be a constant weight binary {0, 1} code of maxi- 
mum possible size having the following parameters: length= 
A, weight= u>, and maximum inner product between any two 
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codewords < k. The size of C cw 
one-dimensional Johnson Bound: 



< 



A- 1 



1 



A 



is upper bounded by the 



A 



The idea is to construct a 2-D OOC whose code arrays are 
partitioned into | C cw | subsets with each subset associated 
to a distinct codeword in C cw . Consider a codeword in C cw 
where the l's in this binary codeword, appear in the ui symbol 
locations Ai, A2, ■ ■ ■ , X w . We associate with this codeword, a 
maximal collection of 2-D code arrays with MCP k which are 
such that only the wavelengths associated to rows Ai, • • • , A„, 
contain a pulse. No pulse is sent along any of the other 
wavelengths. For any of the | C cw | choices of uj wavelengths, 
let us use a 2-D OOC C, 



oppw 



with exactly one pulse per 
wavelength. 

It is easy to see that the composition of these two codes 
forms a 2-D OOC code with parameters (AxT,u,k) of size 

I C C W I • I C ppyj J . 

Since the new code is AM-OPPW, the autocorrelation is 
by Remark [T] 

For crosscorrelation, consider two codewords C a and C& 
from the new code. There are two possibilities: either both 
C a and Cfc have originated from the same codeword (of the 
constant weight binary code), or from different codewords. If 
C a and C correspond to the same codeword in the constant 
weight binary code, then they have the same selection of rows 
and hence, their crosscorrelation is the same as that of the 
corresponding codewords from the 2-D OPPW code, which we 
know is < «. If, on the other hand, C a and Cb have originated 
from two different codewords in the constant weight binary 
code, then the maximum number of rows that can overlap 
in these two codewords is bounded above by ft, and each 
row can contribute a maximum of 1 collisions; hence, the 
crosscorrelation is bounded above by ft. 

Remark 6: One might think of using a Steiner system 
instead of a constant weight binary code. The size of a 
S(k+ 1, u>, A) Steiner system (which is the same as a (k+1)- 
(A, ui, 1) t-design) is 



A 

K+ 1 



UI 
K+ 1 



(24) 



This reduces to 

A / A 

UI 



UJ - 1 



A 



(25) 



Note that every Steiner system is also a constant weight 
binary code and hence satisfies the Johnson Bound for 1-D 
codes with equality. However, Steiner systems do not exist for 
all possible values of k,A,uj. Hence, using constant weight 
binary codes provides a more general construction. 

A. Construction CP1: AM-OPPW, k < uj < A, T is prime 

For the case when T is prime and uj < T, a code 
C ppw of maximum-possible size T K can be constructed using 



construction PI. The overall size of the new 2-D OOC in this 
case is given by 



A- 1 
uj - 1 



A-2 



A - ft 



Proposition 6: In the above construction, the resulting AM- 
OPPW 2-D OOC is asymptotically optimal, if the constant 
weight code is asymptotically optimal (in Johnson Bound). 
For the proof, we refer the reader to Appendix |VII| 
Example 2: In this example, we construct a (7 x 5,3,1) 
AM-OPPW 2-D OOC. We first need to choose a constant 
weight code of length 7, weight 3 and k = 1. We know that 
the Singer construction [2], [62] for OOCs has the desired 
parameters and that its corresponding constant weight code is 
optimal and consists of the following codewords: 

{0,1,3},{1,2,4},{2,3,5},{3,4,6}, 

{0,4, 5}, {1,5, 6}, {0,2, 6} 

Using construction PI for OPPW 2-D OOC, we construct a 
(3 x 5, 3, 1) OPPW 2-D OOC of size 5. Concatenating these 
two codes will result in a (7 x 5, 3, 1) AM-OPPW 2-D OOC 
of size 35, which is optimal by Proposition [3] 



B. Construction CR1: AM-OPPW, uj < A, 

T = p rn + 1 



ft < uj is even, 



In a manner similar to the one described above, we compose 
the OPPW construction Rl with a constant weight binary code. 
The overall size of the new 2-D OOC is given by 



C l< 



c(«0 
T 



1 



A-l 

UJ - 1 



A-2 
uj - 2 



A 



Proposition 7: In the above construction, the resulting AM- 
OPPW 2-D OOC is asymptotically optimal, if the constant 
weight code is asymptotically optimal (in Johnson Bound). 
Proof: Using equation ( |63"j ) of Appendix [VJ we see that 



lim \C\ > T K , 



(26) 



which is the same as the number of codewords in construc- 



tion PI. The rest of the proof follows from Appendix VII 



_L 



v> t 2 t 3 t 4 t 5 t 6 t 7 



Fig. 3. Example of Polarization Rotation Invariant OOC 
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VII. Codes Exploiting Polarization 

A third dimension that can be exploited to construct OOCs 
is the polarization dimension. Light propagates along two 
orthogonal polarization states. Under ideal conditions, these 
two polarization states will be perceived as being orthogonal 
at the receiver [63] despite polarization rotation. 

If we can design codes with good correlation properties that 
are resilient to both timing and polarization ambiguity, we can 
use these codes, which we shall refer to as 3-D OOCs [61], 
[64], to spread the signal in time, wavelength and polarization. 

A 3-D (2xAxT,w,k) OOC C is a family of {0,1} 2xAxT 
arrays of constant weight w. Every pair {A, B} of arrays in 
C is required to satisfy: 



T-l 



2 A 
p=l A=l t=0 



>t r 2 )) < K (27) 



where either A ^ B, t\ ^ or r 2 7^ 0. All the bounds of 



Section III can be generalized to this class of OOCs. 

Polarization Rotation Invariant Codes Construction: It 

is of course possible to construct a 3-D OOC by starting from 
a 1-D OOC and applying the Chinese Remainder Theorem 
(CRT) [65]. When the 3-D code is constructed from a 2-D 
code, however, then the transformation can be seen as a means 
of reducing the required chip rate. 

Let C be a (A x (2T),w,k) 2-D OOC where T is an odd 
integer. For every codeword in C, apply the CRT mapping 
to each wavelength to spread that wavelength in time and 
polarization. The CRT mapping is a one-to-one mapping 
between a sequence and an array that preserves correlation 
values. It follows that by making use of the two orthogonal 
polarization states, we have in effect, reduced the needed chip 
rate by half. 

Fig.[3]shows an example of a 3-D OOC with 3 wavelengths, 
7 time slots, weight 6 and k = 1 constructed using the above 
construction. 

VIII. Sequences for Phase-Encoded OCDMA 

We now focus our attention to designing sequences for 
phase-encoded OCDMA. 

The spectral encoding OCDMA system is harder to imple- 
ment in comparison with direct-sequence OCDMA. That is 
perhaps the reason why most studies on spectral encoding 
systems have an implementation focus. There is not much 
literature on the subject of spreading sequence design with 
the exception of a few results on spectral amplitude encoding. 
Most experimental results reporting on spectral phase-encoded 
OCDMA have assumed synchronous systems and the use of 
Walsh-Hadamard sequences as spreading sequences. As will 
be shown, Walsh-Hadamard Sequences are indeed ideal for 
the synchronous case but quite unsuitable for asynchronous 
systems. Other papers in the literature have proposed the use 
of m-sequences or Gold sequences as spreading sequences but 
do not provide adequate justification for their use. 

We first present a model of an asynchronous phase-encoded 
OCDMA system, and then identify a metric reflective of the 
amount of cross-correlation (other-user interference) in the 
system. Based on this model, we formulate the sequence 



Mode-locked laser 
with K locked 
modes 



K tunable mkroresonators used as 
phase modulators for encoding the 
signal 



OOK 
Modulation 



OOK 




DeModulation 





time 



Shared 
Channel 



Fig. 4. Phase Encoding OCDMA system with Coherent Source. 



design problem. As will be shown, this problem is closely 
related to the PAPR problem in OFDM [66]-[69]. Finally, in 
the next section, generalized bent functions [70] are used to 
construct efficient spreading sequences for an asynchronous 
system. 



A. System Model 

The system that we model in this section is a phase encoding 
OCDMA system with coherent laser source. A diagram of this 
system with a transmitter and receiver is shown in Fig. |4] The 
typical laser sources used for coherent transmission are mode 
locked lasers (MLL). The electrical field of a mode locked 
laser can be written as 



E 



MLL 



(*) 



K-l 

— giwot ^ ' g-ifc(Aw)t 
fc=0 



(28) 



In this equation, K is the number of modes in the mode locked 
laser, and Ao; is the channel spacing between two consecutive 
modes in the mode locked laser. 

The output of the MLL is then passed through a phase 
encoder. In our model, the phase encoder applies different 
phase shifts to different modes of the MLL to spread it. 
Conventionally, the phase masks used in this approach consist 
of only {0,7r} phase shifts. Recently, Stapleton et al. [71], 
[72] showed that by using microdisk resonator technology, 
any phase can be applied to the different modes of the MLL. 
In light of this result, no restriction on the choice of phases is 
considered in this paper. The output of the phase encoder is 
of the form 



E Enc {t) = e^' 1 Y, 



K-l 

iu>ot \ gi(fc(Au;)i+<^s,) 



(29) 



fe=0 



where is the phase shift that the encoder applies to the fc-th 
mode of the MLL. Upon OOK modulation with data bit d 



K-l 



E Tr (t) = de luJot J2 e'f'^'^'.de {0, 1} . (30) 



fe=0 



At the receiver, the phase decoder applies the inverse phase 
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shift — (f>f. to each mode k of the received signal 

K-l 



E 



Dec 



(t) = de luJat 



= de 



lUjQt 



fc=0 
K-l 

£< 

fc=0 



,i(k(Au))t+4> k -<p k ) 



,ife(A(»;)i 



(31) 



When = 0, there is no interference and we are back to 
the single-user case. Hence, we assume S n) = 1 from now. 

Remark 8: In the synchronous r = case, if $( m ) and 
$(«) are Walsh-Hadamard sequences, (i.e., each {exp(i0^)} 
is a sequence in a Walsh-Hadamard sequence family), the two 
summations in ( |33| ) become zero, and there is no interference. 
This is of course clearly not the case when (see Example 



[3l 3.11(1 Fi [si 

which is the original signal in ( |28j ) modulated by data bit d. g ettm 
After the phase decoder, a photo detector is used to detect the 6 m ^ 



intensity of the received signal 



K-l 



\E Dec {t)\ 2 = 



de 



lUjQt 



K-l 

£■ 

fc=0 



= d 



K-l 

£« 



,ik(Auj)t 



B„ m (r) = e 



-^(A^t+W'™'-^)] 



(32) 



fe=0 



and noting that 



«e( e -^ r e nm (r)) |< |6 nm (r)| , 



(34) 



(35) 



we obtain 



d (m) R 2 + |e nm(r )|2 _ 2d (m) K \Q nm ( T )\ < ^ 

< d^X 2 + |9 nm (r)| 2 + 2d<™># |6 nm (r)| . 



(36) 



If we sample this signal at time t = 0, then the received 
signal will be dK 2 and we can retrieve the transmitted data d 
using a threshold detector. 

Remark 7: In this model, no noise source is considered. 
This is because we wish to focus on the effect of multiple 
access interference (MAI). 

When there is more than one user transmitting data, the It follows that minimization of |9„ m (r)| for all r is a 
receiver receives the superposition of the signals. Assume reasonable criterion for signal design. 

that users m and n are transmitting data simultaneously Remark 9: The above generalizes in a straightforward fash- 
and asynchronously. Each user uses its own phase encoder ion to the case of more than 2 users. 

= {(pQ,<f>{, ■ ■ ■ ,4> e K _ 1 } where i G {m,n}. Let the time 
difference between user m and n be denoted by t (t = in 
a synchronous system). The received signal is given by 



B. Connection with PAPR problem 



K-l 



E ( ™\t) + E${t + r) = d^e luJot £ e W*)t+^) + 

K-l 

d (n) e w (t+r) e i(k( 



Our objective is thus the design of sequences of length K 
such that: 



k=0 



t(A U )({+T)+^"') 



max 



K-l 

£< 

fc=0 



- l [fc(A W )r+(0("»-0< m »)] 



(37) 



fc=0 



is minimized for every sequence pair {</>!}> {^i }■ Equiv- 
The signal at the output of the phase decoder tuned to user alently, we seek to minimize 



m takes on the form 

K-l 

d (m) e iu t ^ gj(^(Aw)t) 



max 

tG[0,1) 



K-l 



7i — ik27TT—i((j) 



(38) 



fc=0 



d (n) elWo (t+r) 



o i(fc(A W )(t+r) + (^<"»-0< m) )) 



The design of sequences with minimum PAPR (peak to 
average power ratio) crops up in conjunction with signal 
design for OFDM systems [66]-[69], [73]. Since designing for 
low PAPR is hard, the common design approach is to design 
The output of the photo detector of this receiver wiU be the for low pMEpR (peak . to . mean envelope power ratio), which 



fe=0 



is more tractable. The PMEPR problem (see [69]) is one of 



square of the magnitude of the above expression. As can 

be seen, there is multiple access interference (MAI) at the designing seque nces {a k } that minimize : 
receiver output. Each transmitter-receiver pair is assumed to 
operate synchronously, and consequently, the receiver samples 
its output at time t = to get 



max 

rG[0,l) 



K-l 

£ 

k=0 



a k e 



^ — ik2nr 



(39) 



A = 



K-l 



d {m) K + d {n) e iuoT e4(fe 

k=0 
K-l 

= d^K 2 + d^ 

( K-l 

2d (m) S n) KV\z e~ luJ « T ^ 



l (k(&u)T+(4> i k rl) -<t> k m) )) 



i(fe(Ao;)r+(0' t " ) -0< m) )) 



As can be seen, in our problem, we are interested in phase 
sequences $ (m) and $( n) such that exp(-i(<I> (n) - $( m ))) is 
a sequence with good PMEPR. 

The results in [69] as applied to the present situation are 
stated below. Let 



■ (33) 



fc=0 



Mj = max 

a j=0,- ,K-1 



K-l 



e- lk2 <i)e~ 



*L n) -^ m) )] 



fc=0 



(40) 
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Time, % 



Fig. 6. An Example of application of Proposition |12| with K = 7. 



Fig. 5. An Example of Walsh-Hadamard Sequences with K = 8. 



and 

M„ (if) = max 

T6[0,l) 

From [69], we know 



K-l 

fc=0 



K<t>. 



(»)_,<»>)] 



Theorem 8: M { d K) > 
Theorem 9: For ii" > 3: 



K 



In A' + 0.603 



in which 



1 

6K 



< 



< 



, M C (F K ) 

max r 

F K (t)\M d (F K ) 

]nK + 1.132- 



3 



(41) 



(42) 



and g a complex-valued function defined on Z™. The Fourier 
transform of g is then defined to be the function G given by: 



G(A) 



£ g{x)w-**, A e Z£ 



(44) 



Definition 2: [70] A function /, / : Z™ — > Z 9 is said to 
be bent if all the Fourier transform coefficients of w* have 
unit magnitude. 

Theorem 10: [70] Every affine or linear translate of a bent 
function is also bent. 

Theorem 11: [70] Let q be odd. Then the function / over 
Z q defined by: 

k 2 



f(k) 



ck all k S Z a 



(45) 



is bent for all c £ Z q . 



K-l 



K-l 



F K (t) = a k e 2mkt such that ^ |a fc | 2 = K (43) 
and 



fc=0 



fc=0 



M d (F K ) 
M C (F K ) 



max 

j=0,...,K- 



The implication of Theorem [9] is that, if we design se- 
quences with good asynchronous properties for all the -4= 
samples of r, it is guaranteed that the same sequence has 
good asynchronous properties for all values of r. 



IX. Sequence Constructions based on Bent 
Functions 

In this section we use generalized bent functions to design 
sequences with good asynchronous properties. Some prelim- 
inaries on generalized bent functions that we will use are 
introduced in the first subsection. 



A. Generalized Bent Functions 

Definition 1: [70] Let Z™ denote the set of m-tuples with 
elements drawn from the set of integers modulo q, w = e*("«") 



B. New Construction 1 
Proposition 12: Let 

$ = {*W|* W = {^^»-",^- 1 }^€{l,2,..- ,L}} 

be a family of phase sequences such that the difference 
sequence is associated to a bent function, i.e., 

2tt 



_ $( m ) 



K 



if (0),f (!),-■■ ,f(K-l)),n^m 



where f(x) is a bent function over Z#. Then max|0„ m (r)| 
is as small as it can possibly be over multiples r of K 2 ^^ , 
and thus these phase sequences are suitable for use in asyn- 
chronous phase-encoded OCDMA systems. 



For the proof, we refer the reader to Appendix VIII 



Proposition 13: The following phase sequences have the 
minimum possible value of max |6 n m( r )l property to be used 
for asynchronous phase encoding OCDMA systems with K 
modes, where if is an odd prime: 



j(m) 



(/c 3 + a m k 2 + b rn k + c m ) 



2tt 



G Z K ; m ^ n : a m ^ a„; Ka. prime > 2. 



For the proof, we refer the reader to Appendix IX 
Example 3: Fig. [5] shows the application of Walsh- 
Hadamard sequences for asynchronous systems. In this graph, 

K = 8, u = f, Au = ^, $W = (0,0,0,0,0,0,0,0) and 
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$("■) = (7T, 7T, 7T, 7T, 0, 0, 0, 0). In this figure, the output of the 
MLL as it is seen after the photo detector is denoted by the 
graph that has a maxima above 60 (this curve is red colored in 
the soft copy). The other curve with a maxima just above 30 
(which is blue colored in the soft copy) shows |0 rim (r)| at 
the output. As can be seen, the system has no interference for 
r = (synchronous case). However, even for small deviations 
from t = 0, |6 ram (T)| 2 increases significantly. Because of the 
high peak of |6 nm (T)| 2 , these phase sequence are not suitable 
for asynchronous transmission. 

Example 4: Fig. [6] shows the application of the construction 
of Proposition fT3| where K = 7, gjq = J , Aui — f^, 
(a m ,b rn ,c rn ) = (2,5,3) and (a„,6„,c„) = (5,4,1). For this 
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70 
60 
50 

¥ 40 

CD 

I 30 
20 
10 




-60 



system $< m ) = ^,0,0,0, ^) and 

$ („) = ^ ^ 4^ 2k ;0) io^ 2* j_ In tMs figure ^ the output 

of the MLL as seen after the photo detector is denoted by the 
curve that has its maxima around 50 (red colored in the soft 
copy). The other curve with maxima around 15 (blue colored 



-40 -20 

Time t- 



Fig. 7. An Example of Application of Proposition 1 1 5 1 with q = 9. 



Corollary 16: Let q = H T j , _ 1 pj r i such that p m i n 



in the soft copy) shows |6 nm (r)| at the output. Here, the min{p :) }" =1 , where pj is a prime number. Then the maximum 



circles are samples of |0 nm (r)| for r 



2nj 



As can be size is p r , 



1. 



Proof: Suppose maximum size > p ri 



seen, all these values are equal to 7. It can be observed that Then c$ and cj will exist such that 



|@nm( T )| is l° w for all values of r and the phase sequences 
are thus applicable for asynchronous transmission. 

C. New Construction 2 

The second construction is based on the following family 
of bent functions. 

Theorem 14: [74] Let q be an integer which is neither the 
product of distinct primes nor equal to 2 modulo 4. Then the 
function /(■) over Z g , defined by 



f(k + 1) = f(k) + a k , Vfc, a k G Z„ /(0) G Z q 
is bent if the integers a k satisfy the dual conditions: 



K-l 



a k = (mod s), 



fe=0 



a k + C\ns (mod q), Vfc £ Z s ,Vn G Z 



g/s 



(46) 



(47) 



(48) 



Cj = Cj (mod p mi „). 

Hence, we can say that p m i n divides Ci — Cj. Hence, Ci — Cj 
is not relatively prime to q, which is a contradiction. ■ 
Corollary 17: Let q — p r , where p is a prime and r > 2. 
Then the maximum size is p — 1. 

Proof: Assume that the maximum size > p. Then Ci 
and cj will exist such that 

Ci = Cj (mod p) . 

Hence, we can say that p divides Ci — Cj. Hence, c; — Cj is 
not relatively prime to q, which is a contradiction. ■ 
Example 5: Fig. [7] shows the application of the construction 
of Proposition [15] where p = 3, r = 2, s = 2, q = p r = 9, 



10' 



(a ,ai,a 2 ) = (1,3,5) and (60,61,62) = 



where ci is any integer relatively prime to q and ,s is any 
integer greater than one that has the same parity as q and 
whose square divides q. 

Let g(k) be a second bent function over Z g such that 

g(k+l) = g(k) + b k , Vfc, b k eZ q ,g(0) eZ q (49) 

and 

K-l 

h = (mod s), (50) 



(3,5,7). For this system, /(0) = 1 and g(0) = 1. In this 
figure, the output of the MLL as it is seen after the photo 
detector is denoted by the graph that has its maxima around 
80 (colored red in the soft copy). The dotted curve with its 
maxima around 20 (colored blue in the soft copy) shows 
|@nm(7")| 2 at the output. Here, the small circles around 10 are 



2nj 
9(Au) ' 



As can be seen, all 

2 



fe=0 



bk+ns — h + c 2 ns (mod q), Vfc S Z s , Vn e Z q / S , (51) 

where c 2 is any integer relatively prime to q. 

Proposition 15: Let h(k) = /(fc) — g(k) be a function over 
Z 9 . Then h(k) is a bent function in Z q if c\ — c 2 is relatively 
prime to q. 

For the proof, we refer the reader to Appendix [X] 



the samples of |0 m „(r)| for r 

these values are equal to 9. It can be observed that |0 m „(r)| 
is low for all values of r and hence these phase sequences are 
good for asynchronous transmission. 



X. Conclusion 

In this paper, we presented 9 families of 2-D OOCs. One 
of these families is optimal and the rest are asymptotically 
optimal with respected to the Johnson bound or with respect 
to the new bounds proposed in this paper. A novelty of 
our constructions is the large size. This was achieved by 
constructing optimal families for large values of the MCP since 
the optimal family size increases exponentially in the MCP. 
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Appendix I 
Nonbinary Johnson Bound 

(Proof of Proposition |2| 

The proof is along the lines of the proof of the Johnson 
bound in the case of binary codes. 

Assume a constant weight (A, ui, k) code C of size 
At(A,lu, k). If we arrange all the codewords along the rows 
of a matrix, the total weight of the matrix is ujAt(A, ui, n), 
and thus each non-zero symbol is repeated on an average 
ujA t (a,u,k) t j mes T nuSi there is a symbol a which occurs 

at least " a t times, and this symbol is repeated on an 



average 



AT 



times in each column. It follows that there 

uAt (A,u;,k;) 
AT 



exists a column c with at least — T ^ T : ' ; occurrences of the 
symbol a. However, the number of occurrences of a in column 
c cannot exceed At (A — 1, u> — 1, k — 1) . This is because if we 
select all the rows containing a in column c, and then delete 
this column c from all these rows, we will obtain a constant- 
weight code of length A — 1 and weight uj — 1 with Hamming 
correlation < k — 1. 

It must therefore be that: 

uMt(A, ui, k) 



AT 

A t {A,u,k) < 



< At (A - 1,oj- 1,k- 1) 



AT 



A t (A-1,(J-1,k-1) 



(52) 



By repeating this procedure recursively n times, we arrive 



at: 



TA 

UJ 



A t (A,uj, k) < 
T(A-l) 



(53) 



u> - 1 



r(A - k + 1) 

Ul — K + 1 



At(A — K, Ul — K, 0) 



The proof is then completed by noting that 

(A - k)T 



A T (A - k,uj - k, 0) < 



UJ — K 



(54) 



To see this, let us arrange all the codewords of a (A — 
k,uj — k, 0) constant weight code C along the rows of a matrix. 
To satisfy the constraint of zero Hamming correlation, each 
alphabet can occur only once in each column. Thus, there 
are at most T non-zero entries in each column. Since there 
are A — n different columns, the entire matrix can have at 
most (A — k)T non-zero entries. On the other hand, each row 
has exactly uj — k non-zero entries, so that there are at most 

^z^ j rows in the matrix and the assertion is proved. 

Appendix II 
Bound on AM-OPPW Code Size 

(Proof of Proposition [3} 

Let C be an AM-OPPW 2-D OOC of size $(A X T,u,n). 
Create a code C' that consists of all T columnar cyclic shifts 
of each code in C. This code is of size <1>'(A x T,oj,k) = 
T$(A x T,oj,k). 

By identifying each row of a code matrix in C having a 1 in 
the t-th column with the symbol t belonging to {1, 2, . . . , T}, 
and a blank row with the symbol 0, we obtain from the 2-D 



OOC, a 1-D constant weight code over the alphabet Zt+i- 
This 1-D constant weight code has parameters (A, gj,k) and 
is of size <1>'(A x T,uj, k) over an alphabet of size T + 1. 
It follows from our bound above in Proposition [2] that 



$'(A x T,u,k) < 



$(AxT,u,k) < 



TA 

UJ 



T(A-l) 



T(A - k) 



1 


TA 


T(A - 1) 


T(A - k) 


f 


UI 


OJ-1 


OJ — K 



< 



T(A-l) 



T(A - k) 



Appendix III 
Asymptotic Optimality of construction P2 

The number of codewords in this code are 



d|(A-l) 

Consider the largest term in the summation - this corre- 
sponds to d = 1 and evaluates to 

A K+1 _ I 



T 

Hence, the total number of codewords \C\ is given by 

A K+1 - 1 1 



\C\ = 



T 



T 



j: (A\^-i)Kd) 



d|(A-l), d>l 

For A and T tending to infinity, we get 

A K+1 

lim ICI > . 

A,T^oo 1 1 ~ T 



(55) 



Since to = T for this construction, the two-dimensional 
Johnson bound given by Theorem [T] specializes to 



AT - 1 
T — 1 



A (AT- 1 / AT — 



$(Axr,w,)t) < 

< 

- T \ T —1 
For A and T tending to infinity, we get 



AT -k 
T-k 



V T-K 



lim $(A x T, us, k) 

A,T^oo 



< 



A /AT 
T \ T~ 

A K.+ 1 

T ' 



AT 



(56) 



From ( |55] > and ((56]), we can see that construction P2 is 
asymptotically optimal. The proof for asymptotic optimality 
of constructions P4 and P5 is along similar lines, since the 
expression for the total number of codewords is similar. 
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Appendix IV 
Asymptotic Optimality of construction P3 

The number of codewords |C| in this code are 

(T+ 1) K+1 - 1 



> 



T 

T ' 



When A and T tend to infinity, we get 

lim \C\ > 



(57) 



Since this construction is AM-OPPW, we use the bound 
given by Proposition [3] to get 



$(A x T,w,k) < 



A 



T(A-l) 

to - 1 



T(A - k) 



< 



A /T(A-l) / T(A - k) 

W\LJ— 1 \ — K 

A / T(A - 1) / T(A - k) 



A — K \ A — K — 1 \A — K — K 
since lo = A — n. When A and T tend to infinity, we get 



A / TA 

lim $AxT,(j,k < - — 

A.T^oo " A \ A 



TA\ 



(58) 



Comparing ( |57| i and ( |58] l, we can see that this construction 
is asymptotically optimal. 

Appendix V 
Asymptotic Optimality of construction Rl 

The number of codewords |C| in this code is 



where 



c(k') = 



> 2k' + 1 _ <f^_ 



k' = 1,2,3,4,5,6 
k' > 7. 



We consider the two cases separately: n' < 6 and k' > 7. 
For k' < 6, we have 



\C\ 



> 



£M +1 
T 

T 

q 2K ' +1 - q 
T 

(T-1) 2K ' +1 



(59) 



T + 1 



T 

In the limit A, T tending to infinity, we get 



lim \C\ > T 



2n' 



(60) 



For k' > 7, we have 



\C\ = ^+1 



> 



> 



T 

„2«'+l _ g 2 "' -6 
^ 7 



(61) 



T 



(T-l) 



2«'+i _ (r-i) 2 

7 



In the limit A, T tending to infinity, we get 



lim \C\ > T 2k ' 

A.T^oo 

= T K . (62) 
We can see that ([60]) and ([62]) are the same. Hence, 

lim \C\>T K (63) 

for all values of k', i.e., for all even values of k. 

The two-dimensional Johnson bound given by Theorem [T] 

is 



$(AxT,w,k) < 



A 


AT 


- 1 


LO 


LO - 


- 1 




AT - 


- 1 


= ( 


LO — 


1 




AT - 


- 1 



AT -k 



LO 


— K 


AT 


— K 


LO - 


- K 


AT 


— K 


A- 


- K 



A V A- 1 

since ui = A for this construction. In the limit A, T tending 
to infinity, we get 



lim $(A x T.w.k) < T K . 

A.T^oo 



(64) 

By comparing ( [63) and ( |64[ i, we can see that this construc- 
tion is optimal. 

Appendix VI 
Asymptotic Optimality of construction R2 

The number of codewords |C| in this code is 

M(d, T) 



\C\ 



T 



h(x) 



(i{h{x))N{d- deg(h(x)),T). 



Note that N(-) is an increasing function in its first argument, 
which gets maximized when the degree of h(x) is 0. Moreover, 
for h(x) = 1, jj,(h(x)) achieves its maximum value of 1. Since 
the summation above is greater than the largest term of the 
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summation (which corresponds to h(x) = 1), we conclude 
that 



\C\ > -N(d 7 T). 

Since we are interested in proving the asymptotic optimality, 
it is sufficient to work with the highest order term. Substituting 
the value of N(-), the above sum simplifies to 



T(q-1Y 



\C\ > 



From the definition of u{\) in ( p2) , we notice that its 
maximum term corresponds to I = 1. This is so because 
y(-) decreases exponentially with I. Moreover, attains its 
maximum value of 1 for 1=1. Substituting the value of (the 
largest term in the summation of) u(l) from |22) into the 
above summation, we get 



\C\ > 



2/(1) 



T(q-1) 



(65) 



A- 1 

UJ- 1 



A-2 
u-2 



A- 



K 









( A 


EH 






1 


V UJ 





UJ — K 
A - K 



> 



- 1 



lim 



A-l 

UJ - 1 



UJ — K 

A-2 
uj - 2 
A K+1 



-1-1-1 



A - K 

UJ — K 



> 



uj(uj — 1) • • • (u) — k) 
=> lim | C cw | T K > -. 

A^oo UJ[UJ — 1) ■ ■ ■ (U) — K) 



On the other hand, from the bound in Proposition [3j 



$(A x T,u,k) < 



T(A - 1) 
uj - 1 



T(A - k) 



< 



AT(A-l) T(A-k) T k A(A-1)---(A-k) 



w u; — 1 



w(w — 1) • • • (uj — k) 



We now compute an approximation for y(l) and then 
substitute it back in the above inequality. We do this by 
substituting I = 1 in ( |22"] > to get the following approximation 
for y(l) 



2 g { g rf- ei +i { q d -^+ c+i - {q d - ei+l - iy 



=> lim *(A x T,uj,k) < - 

A^oo ujyuj — 1) ■ ■ ■ (uj — k) 

which establishes what we need to prove. 



Appendix VIII 
Proof of Proposition[T21 



ei=0 



31=0 



K-l 

E 

fe=0 



( 2ttj 



K-l 



\K(Auj) 



e -*^j'fc e -*1r/(fc) 



fc=0 



E 

fc=0 



= F(j), 



Each summation in this expression is a geometric progres- 
sion. Summing the series and approximating (by neglecting 
the lower-order terms in the summation), we get 

+ q 2d+2 

Substituting this value of y(l) into ( |65] >, we get 



where w = e l £ . Since the function f(x) is a bent function: 
1 



so for all j, 



F(j) = 1 => F(j) = VK 



pi Z' \ 



is the minimum value 



a, ( \ 
^~>nm\ k(Au>) > 



= a/A 7 , and by Theorem |sj this 
^U— ^ can take. 



By Theorem [9] this sequence is good for all values of r. 



\C\ > 



3q 



2d+2 



T(q-l) 
3(A - 1) K+2 



TA 



> 



A «+i 
T ' 



(66) 



Comparing this with ( p6[ l, we see that construction R2 is 
asymptotically optimal. 



Appendix VII 
Asymptotic Optimality of Construction CP1 



(Proof of Proposition |6| 



Appendix IX 
Proof of Proposition[T31 

For any two different sets of phases $( Il ),<I)( rn ) from the 
given set of phase sequences, we have 

2tt 



and 



6 M 



(fc 3 + a m k 2 + b m k + c m ) 
(k 3 + a n k 2 + b n k + c„) 



K 
2tt 



So 



((0„ - a m )fc 2 + (6„ - & m )fc + (C„ - Cm)) 



2tt 
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Since we assume a n ^ a m , it implies that 

(a n — a m )k 2 + (b n — b m )k + (c n — c m ) 



From equations |47) and (|50j), we know that both these sums 
are equal to 0, and so is the difference. 

We next need to show that p k satisfies the second condition: 



is a bent function by Theorems 10 and 11 So the given set 



of sequences satisfies the conditions of Proposition 12 and is 
good to be used with asynchronous phase encoding sequences. 

For any two different set of phases $(*),<!>(■?) from the above 
set of phases, we have: <j>^ = (fc 3 + a,;fc 2 + 6;fc + Cj)^- and 



Pk+ns Qk+ns b k j rns . 



(69) 



t>f> = (fc 3 + ajfc 2 + b j k + c J ) 2 ^. Then 



2tt 



4>f ~ <Pi j) = (K - %")fc 2 + (6* ~ b 3 )k + (a - Cj)) -—. (67) 
Since we assume a, ^ aj, it implies that (<z,j — aj)k 2 + (hi — 



bj)k + (c, — cj) is a bent function by Theorems 
Hence, 



10 



and 



11 



K-l 



1 ST „(«*-. 



j )k 2 +{bi-b j )k+(c i -c j ) \k 



fc=0 



= e i (^),A = 0,l,...,i<r-l 



K-l 



e i[(«i-ai)fe 2 +(6i-^)fe+(c«-c ;j )](^) e i(^)fc 



k=0 



= VK 



fc=0 



K. (68) 



So |6 nm (r)| for all r = ^ is equal to \/K. Using 
Theorem |i] we have met M^ K \ In addition, using Theorem 
[9] we can conclude that |0 nm (r)| is almost equal to \[K for 
all values of r. So the sequences described above are good 
asynchronous phase encoding OCDMA sequences. 

Appendix X 
Proof of Proposition[T31 

We know that 

h(k) = f(k)-g(k) 
=>h(k+l) = /(fc+l)-. 9 (fc + l). 

Using equations ( |46] > and (^9j, we get 

h(k + l) = (f(k)+a k )~(g(k) + b k ) 
= (f(k)-g(k)) + (a k ~b k ). 

Hence, we can say that 

h(k+l) = h(k)+ Pk , 

where p k = a k -b k . 

For h(k) to be bent, p k should satisfy the dual conditions. 
Consider 



K-l 



K-l 



k =0 k=0 

K-l K-l 



k=0 



fe=0 



Using equations ( |48| ) and ( pi) , the above equation can be 
rewritten as 

Pk+ns = (an, + cms(mod q)) - (b k + c 2 ns(mod q)) 

= (a k -b k ) + (ci - c 2 ) ns(mod q) 

= p k + (ci — c 2 ) ns(mod q). 

Hence, p k satisfies the dual conditions provided that c\ c 2 
is relatively prime to q. Hence, h(k) is a bent function. Hence, 
the given set of sequences satisfies the condition of Proposition 
12 and is good to be used with asynchronous phase-encoded 
OCDMA. 
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